arXiv: 1501.06509vl [astro-ph.HE] 26 Jan 2015 


Mon. Not. R. Astron. Soc. 000, ^{5] (2014) Printed 27 January 2015 (MN style file v2.2) 


Modeling astrophysical outflows via the unifled 
Dynamo-Reverse Dynamo mechanism 

Manasvi Lingam^* and Swadesh M. Mahajan^f 

^Institute for Fusion Studies, The University of Texas, Austin, TX 78712, USA 


27 January 2015 


ABSTRACT 

The unified Dynamo-Reverse Dynamo (Dy-RDy) mechanism, capable of simultane¬ 
ously generating large scale outflows and magnetic fields from an ambient microscopic 
reservoir, is explored in a broad astrophysical context. The Dy-RDy mechanism is 
derived via Hall magnetohydrodynamics, which unifies the evolution of magnetic field 
and fluid vorticity. It also introduces an intrinsic length scale, the ion skin depth, al¬ 
lowing for the proper normalization and categorization of microscopic and macroscopic 
scales. The large scale Alfven Mach number A4a, defining the relative “abundance” 
of the flow field to the magnetic field is shown to be tied to a microscopic scale 
length that reflects the characteristics of the ambient short scale reservoir. The dy¬ 
namo (Dy), preferentially producing the large scale magnetic field, is the dominant 
mode when the ambient turbulence is mostly kinetic, while the outflow producing 
reverse dynamo (RDy) is the principal manifestation of a magnetically dominated 
turbulent reservoir. It is conjectured that an efficient RDy may be the source of many 
observed astrophysical outflows that have Ma^ 
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1 INTRODUCTION 

Astrophysical jets are universal: they appear in myriad 
forms, exhibiting significant variations in luminosity, colli- 
mation and velocity. Relativistic jets have been observed 
in the context of Active Galactic Nuclei (AGNs), micro¬ 
quasars and Gamma Ray Bursts (GRBs). Smaller-scale non- 
relativistic versions have been observed in a much wider 
array of systems such as low-mass Young Stellar Objects 
(YSOs), massive X-ray binary systems and Protoplanetary 
Nebulae (PPNe). Although the word “jets” suggests a highly 
collimated and narrow outflow, diffuse emanations, often 
termed “outflows”, have also been detected in stars in their 
final stages of evolution. In addition, our Sun exhibits jets 
and flares, albeit ones that are much smaller than the one 
mentioned above. In this Letter, we shall refer to jets and 
their (relatively) diffuse cousins as “outflows”. A comprehen- 
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Given the diversity of their occurrence and the 
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vastly different parameters that these outflows ex¬ 
hibit, the quest for a broad underlying mechanism 
that can power these outflows has proven elusive 
although several seminal works in the 1970s and 
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netic field, in their attempts to explain the origin of these 
outflows. 

Most of the preceding analyses used the ideal magneto¬ 
hydrodynamic (MHD) model, the simplest and highly ex¬ 
plored plasma fluid model. However, in the last decade, 
the importance of the Hall effect in astrophysical sys¬ 
tems has been recognized; the Hall effect has been shown 
to play a crucial role in sever al scenarios, such as mag¬ 
netic braking and star formation ( Kras nopolskv. Li fc Shand 
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along with the realization that MHD may be an inad- 
equate framework for generating the observed outflows 
ifBisnovatvi-Kogan fc Lovelacell2001^ . 

In order to explore the simultaneous dynamics of mag¬ 
netic fields and outflows, one must resort to models more 
encompassing than MHD. As an illustration, we will ana¬ 
lyze incompressible Hall magnetohydrodynamics (HMHD), 
the simplest “beyond MHD” model that captures one cru¬ 
cial two-fluid effect: the Hall current, expressing differential 
electron-ion motion, introduces a fundamental scale length, 
the ion skin depth (ideal MHD is scale-free) that provides 
a convenient fiduciary length in terms of which one can de¬ 
fine ‘large scale’ and ‘small scale’ quantities. HMHD assigns 
assigns equal weights to, and thereby unifies, the magnetic 
fields and flows in a natural manner, as demonstrated in 
Section 0 Qualitatively, this can be understood as a conse¬ 
quence of one of the HMHD equations involving a simultane¬ 
ous evolution of the magnetic field B and the fluid vorticity 
V X V, resulting in the inexorable linking of these two quan¬ 
tities. 

In this Letter, we explore what was termed as the 
‘reverse dynamo’ m echanism; derived and developed in 
iMahaian et al.lll2005h . The mechanism employs HMHD, and 
exhibits potential universality in generating non-relativistic 
outflows. The original moniker ‘reverse dynamo’ consti¬ 
tutes an inadequate (and incomplete) description of the 
more encompassing process. Though dyn amo based ap - 
proaches have been around since the 1970s (lLovelacelll976ll . 
with se veral important developments in the subsequent 
decad es (lBhnfTachariee_^£jImneir|[l98^; Brandenburget^^] 
19951: Iv ishniac fc Choi l200ll: iB lackman fc Fieldl 200d; 


Brand«nbuTgJ_ 2005L Brandenburg fc SubramaniarJ 20051 : 


Ebrahimi fc Bhattacharieell20l4 l. the Dv-RDv approach is 
conceptually unique as it incorporates the following features: 


• The HMHD is the governing model from the outset (the 
Hall term is not appended to an essentially MHD calcu¬ 
lation). Consequently the magnetic field and the flows are 
treated on an entirely equal footing. The “dynamo” desig¬ 
nates the creation of long-scale magnetic fields from a short 
scale (turbulent) velocity field, and the “reverse dynamo” 
denotes the complementary/opposite phenomenon of gener¬ 
ating long-scale flows from a short scale (turbulent) mag¬ 
netic field. The HMHD based theory is a unified Dy-RDy 
theory; both flows and magnetic fields emerge simultane¬ 
ously from a given source of short scale energy be it kinetic 
or magnetic. 

• Unlike MHD, the HMHD has an intrinsic length scale. 
Thus it becomes possible to characterize long and short 
scales in a well-defined way. This wifi prove to be crucial 
in estimating the length scale relevant to, for example, the 
generated outflows. 


2 DY-RDY: A SYNOPSIS 

We begi n with a short r ecapit ulation of the formalism de¬ 
rived in IMahaian et ahl (1200 Sl l. Assuming a constant den¬ 
sity (purely for analytic simplicity) and equal temperatures 
for the electrons and ions {p = Pi + Ps ‘2nT), the HMHD 
yields the following evolution equations for the flow and the 


magnetic field: 

^ = V X [(v - aoV X b) X b], 

^ = vx(Vxv) + (Vxb)xb-v(^p+y^. (1) 

wherein ao = ^io/Ro, with Aio representing the ion skin 
depth, which serves as the intrinsic length scale character¬ 
istic of HMHD. Rq represents a characteristic scale length 
to be specified shortly hereafter. In m, the magnetic field 
is normalized to an arbitrary representative magnetic field 
Bo, the velocity to the corresponding Alfven velocity Vao, 
the length scale to Ro and the timescale to Ro/Va- Subse¬ 
quently, we choose Ro to be the ion skin depth, as we are 
interested in the underlying microscopic physics. Thus, we 
observe that Aio will no longer appear explicitly, i.e, ao = 1. 
With such a choice of units , one may re-express (fTI) as a pair 
of vorticity-like equations (IMahaian &: Yoshidalll998l l: 

^ - V X (v, X D,) = 0, 

where the vorticities flj and the associated velocities Vj are 
Di = B, vi = V — V X B, 

122 = B -b V X V, V 2 = V. 


The structure of the above equations is radically different 
from that of ideal MHD, and is also indicative why B and 

V (V X V to be precise) are generated in tandem. 

Notice that the governing equation for the canonical 
vorticity 122 treats the magnetic field and the vorticity on 
par. There are two conserved helicities in the system: the 
magnetic helicity A ■ B and the canonical helicity 

d^x (A -b v) • (V X V -b B); the latter is essentially the 
ion helicity reflecting the (finite) ion inertia. Unlike in ideal 
MHD, the cross helicity, d^xv-B, is no longer conserved. 

The velocity and magnetic fields are decomposed into 
the equilibrium seed fields (vo and bo) and the fluctuations; 
the latter, in turn, comprising of the macroscopic (U and H) 
and microscopic (v and b) components: 

b = H + b -b bo, 

V = U-bv-bvo. (2) 


It is important to realize that the fluctuations are comprised 
of both macroscopic and microscopic components. The for¬ 
mer yield non-zero expressions upon carrying out a suit¬ 
able ensemble or spatial averaging, whilst linear combina¬ 
tions of the latter yield no contributions. Over a microscopic 
scale, we also assume that the large scale fluctuations do not 
vary significa ntly, i.e. that their derivatives vanish. We re fer 
the reader to Mhfimii. Gomez fc MahaianI (1200 2l . I2003al f9l ; 
IMahaian et al.l 1 2005^ for a detailed discussion of the closure 
scheme, and its accompanying assumptions. 

The equilibrium fields are the energy reservoir for driv¬ 
ing the fluctuations. A comment regarding the length scales 
is in order. The only intrinsic scale is defined by the normal¬ 
izing length, the ion skin depth. Any macroscopic length 
scale of interest is way bigger than Aio, but there is some 
latitude in what could be called the microscopic length. The 
latter (microscopic length) is taken to be any length that is 
on the order of, or smaller than, the ion skin depth Aio (or of 
order unity and less in normalized units). The equilibrium 
fields are small scale (microscopic), have been produced by 
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some microscopic process, undergone saturation, and have 
built up an energy reservoir that will drive the large scale 
(macroscopic and observable) and small scale (microscopic) 
fluctuations. 

Given that one has little information about the equi¬ 
librium fields, one makes the most ‘natural’ choice in the 
context o f Hall MHD; the double B eltrami equilibria, de¬ 
scribed in iMahaian fc Yoshidal (Il998l l. The double Beltrami 
equilibria are the solutions of 


bo 


-I- V X bo = Vo, 


bo -I- V X Vo = dvo. 


( 3 ) 


which, along with the Bernoulli condition 
V ^p-|- =0 constitute the stationary solutions of ©• 

The constants a and d, set by the helicities of the system 
dMahaian et al.l l2005li . define two (dimensionless) inverse 


scale lengths A± = 


(d — a ^) ± y^(d -I- — 4 


. Ap¬ 


propriate tuning of a and d can lead to vastly separated 
roots, A+ = A and A_ = ^; the former (latter) represents 
the microscopic (macroscopic) inverse scale length. Such a 
choice is physically meaningful - the astrophysical systems 
are macroscopic, but their underlying physics may have a 
substantial microscopic component. In fact, we could as¬ 
sume that the ambient fields are entirely microscopic; this 
is again a rational assumption since the underlying reser¬ 
voir that generates observable (macroscopic) phenomena is, 
indeed, presumed to be microscopic. 

With these assum ptions and a great deal of algebra 
l|Mahaian et al.l l2005l 'l. one derives the final equations for 
the macroscopic fluctuations: 


H = -r (V X H) , U = V X (sU - qH) (4) 


where r, s and q, given by 

q = r = A^ (l - Aa"^ - a”^) , 

s = a| [(A + a-i)"-l] , (5) 

are determined by b^, representing the normalized ambient 
magnetic energy, and by the scale lengths set by the mi¬ 
croscopic helicities; the entire macroscopic dynamics is con¬ 
trolled by the ambient microscopic dynamics. By manipu¬ 
lating the equations in ®, we find that 

U=—^H. (6) 

s -I- r 

Thus, by solving the first equation in (|4]), we can fully de- 
te rmine H and U . The details of the solution can be found 
in iMahaian et af] (l2005li . 


3 DY-RDY IN ASTROPHYSICAL SYSTEMS 

Equation[5]is the most potent and revealing statement of the 
Unified dynamo - the simultaneous generation of large scale 
magnetic fields and flows out of a short-scale reservoir of 
kinetic and magnetic energy. The dominant dynamo |H| ';$> 
|U| or the dominant reverse dynamo |U| 3> |H| state are 
simply the two limiting cases of the same mechanism. 

What mix of |H| and |U| will finally emerge depends 
upon the coefficients q, s,and r, which are in turn deter¬ 
mined by a and d. The latter duo are also manifest in the 


scale lengths, namely A~^ and The helicity measures, 
in turn, reflect the initial state of the system, in particu¬ 
lar, the ratio of the ambient kinetic energy to the ambient 
magnetic energy. No matter what their ratio, the flows and 
the magnetic fields are generated together; the HMHD is a 
theory of the unified Dynamo-Reverse Dynamo (Dy-RDy) 
process. 

Simple algebra reveals that for a ~ d <C 1, the short 
scale reservoir is mostly magnetic, vo ~ fflbo <C bo leading 
to strong (compared to magnetic fields) macroscopic flows 
U ~ AH > H. This limit pertains only to those systems 
where the observed flows are super-Alfvenic. The standard 
dynamo scenario pertains to the opposite case: a ~ d ^ 1, 
implying the microscopic relation vo ~ obo bo, and the 
macroscopic relation H U. 

This deterministic relationship between the energy mix 
in the short scale reservoir, and the consequent long scale 
observables is one of the major results of the paper. The the¬ 
ory has provided a probe into the earlier era microphysics; 
measuring the long range fields and flows yields information 
about the source that caused them. 

Let us explore a little further the regime where RDy 
part of Dy-RDy will be dominant. We will examine the rele¬ 
vant relation U ~ AH in physical units. Since H and U are 
normalized by Bq and Vao = Boly/Anniin (with constant 
n), we find that 


A 


U 

H / \/ Anmin 



( 7 ) 


Hence, the inverse (dimensionless) microscopic scale length 
A turns out to be approximately equal to the large scale 
Alfven Mach Number Ma, where U and Va are, respec¬ 
tively, the (dimensional) large scale flow and Alfven veloci¬ 
ties. Equation © captures the very essence of this analysis, 
as it directly ties together the microscopic and macroscopic 
physics. This is also seen in as it relates q/(s -I- r) (mi¬ 
croscopic) with U/H (macroscopic). In our subsequent dis¬ 
cussion, we drop the overbars and note that all quantities 
on the RHS of (Q are large scale and dimensional. 

The reverse dynamo mechanism, favored by the order¬ 
ing a ~ d <C 1 (A 3> 1), will operate exactly in those regimes 
where the observed Ma is much greater than unity. Equiv¬ 
alently, if the observed large scale flows are highly super- 
Alfvenic, then it is quite likely that they are generated via 
the reverse dynamo mechanism. The characteristic Mach 
number Ma of a given outflow is an index of the relative 
efficiency of the RDy and Dy mechanisms, which in turn, is 
a mirror of the constitution of the ambient state - what mix¬ 
ture of magnetic to kinetic turbulence it is endowed with. 

Let us now use 0 to probe the range of densities and 
magnetic fields for which the RDy is likely to be important. 
Choosing A > 10 (which loosely satisfies A 1) and differ¬ 
ent values of U, we can calculate Va, and the corresponding 
values of H and n. The results are presented in Fig. [T] For 
the current non-relativistic treatment, we keep U bounded 
below 10"^ km/s. Each of the curves in Fig. [T]corresponds to 
a fixed value of U, and allows us to probe the ranges of n 
and B for which the reverse dynamo may be operational. In 
the regions lying above and on the curve, RDy will dominate 
and generate large scale outflows, while in the regions below 
the curves, the RDy will be sub-dominant. 

In Table [T] the computed Alfven Mach numbers are 
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Table 1. The Alfven Mach number for astrophysical outflows 


System 

p{glcm?‘) 

H{G) 

U{cml s) 

Ma 

Reference 

GRBs 

10^'’ 


10^*^ 

10® 

Beskin (20101 

Microquasars 

lO’^ 

IQi'J 

IF® 

10® 

Beskin (20101 

Radio pulsars 




10“ 

Beskin (20101 

YSOs 

10 

10® 

10' 

10“ 

Beskin (20101 

PPNe 

1 

lo® 

10' 

10“ 

Huffffins ('2007'): Asensio Ramos et al. ('2014') 


The table illustrates the values of the Alfven Mach number A4a for astrophysical systems where 
jets are observed. In the case of Protoplanetary Nebulae (PPNe), the outflows are modeled as 
originating from the central star, which possesses specificat ions akin to the Sun. An upp er bound 
on the large scale magnetic field in PPNe was obtained inlAs ensio Ramos et ^ ll2014h . and the 
typical jet velocities in PPNe were t aken from iHugginsI ll2007l ). The rest of the parameters are 
presented in Table 1 of iBeskinI ||2010| ~). 


n 



Figure 1. The figure depicts the permitted values of n and H 
for the reverse dynamo mechanism. For all the curves, A = 10 
has been chosen. The black, red, green, blue and orange curves 
correspond to the cases with U = 1,10,10^ , 10® , 10^ km/s re¬ 
spectively. Note that this is a log-log plot. 

displayed for a variety of astrophysical systems. The com¬ 
putations should be accepted with the following caveats: 

• The magnetic fields {H) in the large scale outflows can¬ 
not be easily measured, and hence we have chosen to use the 
magnetic fields present in the ‘engine’ (driving the outflows) 
as a substitute. However, since most of the values of A4 a 
are very high, it is reasonable to suppose that reducing the 
magnetic field by several orders of magnitude may still plant 
the system firmly in the regime where the reverse dynamo 
mechanism determines the outcome. 

• The magnetic fields in both the ‘engine’ and the out¬ 
flows are both very complex, and the Alfven Mach number 
serves only as a simplified criterion for evaluating the va¬ 
lidity of the reverse dynamo mechanism. In addition, it is 
known that the winds in GRBs, microquasars and pulsars 
are relativistic, and a fully consistent treatment must take 
such effects into account. We plan to present a relativistic 
version in a forthcoming paper. 

The Alfven Mach numbers, displayed in Table [U are all 
very high - much greater than unity. The simplest criterion 
that such outflows may originate in a RDy mechanism is 
amply satisfied. It must be borne in mind that this does not 
necessarily imply the operation of the RDy mechanism. Nev¬ 
ertheless, this simple criterion can help us probe the nature 
of the underlying source that drives large scale outflows and 
magnetic fields. An example in contrast is the solar wind, a 
highly sub-Alfvenic outflow; the RDy is not likely to be the 
primary driver. However, we note that there are regions in 


the chro mosphere where the R Dy has been shown to play a 
key role jMahaian et al.l[2002l '). 

Lastly, we note that the HMHD is a theory of unified 
Dy-RDy mechanism - the general theory is just as valid in 
scenarios where Ma ^ 1- Even in this sub-Alfvenic regime, 
the flows, along with the dominant magnetic fields, will con¬ 
tinue to amplify as long as there is turbulent energy to drive 
them. Both fields grow at the rate determined by the dis¬ 
persion relation, = — |r| k ilMahaian et al.l 120051 ) with r 
given by ®. One expects the outflow level to be boosted 
over time, albeit in a sub-Alfvenic setting, until it is com¬ 
patible with observations. 


4 DISCUSSION AND CONCLUSION 

In this Letter, we have explored the unified Dynamo-Reverse 
Dynamo mechanism within an HMHD model. The RDy 
(Dy), operates preferentially when the magnetic component 
(the kinetic component) accounts for the bulk of the short 
scale energy reservoir. The end product of dominant RDy 
(Dy) is a super Alfvenic (sub Alfvenic) outflow, where the 
Alfven Mach number obeys Ma f {Ma <C 1). Conse¬ 
quently, if the observed large scale outflows were to satisfy 
Ma 1 {Ma "C 1), they may indicate the operation of a 
reverse dynamo (dynamo) mechanism. A knowledge of Ma 
as per this model can serve as a probe of the “content” of 
the underlying ambient reservoir. 

The unified Dy-RDy mechanism is characterized by a 
striking reciprocity between the micro and the macro scales: 
the ambient (microscopic) Alfven Mach number, denoted by 
Mfj,A, precisely equals {Ma)~^ - the inverse Alfven Mach 
number of the large scale (macroscopic) outflow. 

We have also found that the reverse dynamo mechanism 
may operate in a wide variety of astrophysical objects whose 
ambient density n and the large scale magnetic fields H lie in 
the ranges delineated in Fig.[T]and Table[T] We do not claim 
that the RDy mechanism is the sole source of super-Alfvenic 
outflows. Instead, we emphasize the importance of the large 
scale Alfven Mach number Af a as a means of gauging the 
composition of the ambient magnetic and kinetic energies of 
the reservoir. However, our analysis is not exhaustive; some 
of the inherent limitations were outlined earlier. 

We must stress that the Hall term, which introduces an 
intrinsic micro-scale, and thereby opens up the possibility of 
physics at disparate macros and micro scales, is an essential 
enabler of the unified Dy-RDy theory. This step, however. 

































represents only the tip of the iceberg as Hall MHD, itself, 
is the simplest of the extended MHD models. A natural ex¬ 
tension of the formalism presented in this paper involves 
the construction of an extended MHD and/or two-fluid uni¬ 
fied Dy-RDy that can capture the non-ideal MHD effects 
in a more direct manner. To unders tand relativistic flo ws, 
we must resort to relativistic MHD dLichnerowicz W6^}or 
the coupled relativistic magnetofluid model jMahaianI 2003ll . 
For incompressible Hall MHD, we hypothesize that the re¬ 
sults are likely to remain similar under v —>■ yv in the rel¬ 
ativistic regime. Lastly, we note that the assumption of ho¬ 
mogeneity can also be relaxed without difficulty, and the en¬ 
suing model is still easily solvable via computational means. 

The unified Dy-RDy mechanism presented herein 
(and/or one of its several suggested extensions) constitutes a 
very viable candidate, as well as a strong indicator, for gen¬ 
erating large scale outflows in a wide range of astrophysical 
systems. A more detailed picture of the model’s strengths 
and limitations are likely to emerge only via further numer¬ 
ical simulations. 
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